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Abstract 

Using the new scalar and vector degrees of freedom derived from the non-linear gauge condition 
(d ■ D){d ■ A) = 0, we show that the effective classical dynamics of the vector fields (identified as 
"gluons") in the stochastic vacuum defined by the scalars result in the vector fields acquiring a 
range of possible masses and losing their self-interactions. From this range of masses, we derive the 
mass gap in pure Yang- Mills theory. Finally, we comment on the gauge-invariance of the result. 
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I. INTRODUCTION 



The non-linear gauge condition 

(d ■ D) a \d ■ A b ) = (D ■ d) ab {d ■ A b ) = h(D ■ d) + {d • D)] ab (d ■ A b ) = 0, (1) 

was proposed by the author— as a generalization to the Coulomb gauge because the 
non-linear sector of (1), i.e., those configurations that satisfy d ■ A a = f a (x) ^ 0, cannot be 
gauge transformed to the Coulomb gauged. Subsequent papers established that: 

(a) The gauge condition interpolates between the Coulomb gauge (with transverse fields 
in the short-distance regime) and the quadratic condition (with the new scalar and vector 
fields) in the large-distance regime 3 -. 

(b) The gauge condition can no longer be extended to that with higher orders of the 
Fadeev-Popov operator, i.e., (d ■ D) n (d ■ A) = for n = 2,... is not consistent 3 . 

(c) The effective action of the scalar degrees of freedom is infinitely non-linear and 
a stochastic treatment of a class of classical configurations, those that are spherically 
symmetric in 4D Euclidean space-time, leads to the linear potential and the area law 
behaviour of the Wilson loorA 

(d) The full quantum treatment of the scalar degree of freedom shows that its dynamics 
is equivalent to that of a 2D 0(1,3) non-linear sigma model^. 

In this paper, the effective classical dynamics of the vector field, which we will identify 
as the "gluons", in the stochastic background defined by the spherically symmetric scalar 
will be considered. We will show that the "gluons" will acquire a mass and become non 
self-interacting. 

This paper is then arranged as follows: In Section II, we give the basic equations that 
govern the new degrees of freedom derived from the non-linear gauge. Section III focuses 
on the effective dynamics of the "gluons" and derives an expression for the masses and the 
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mass gap in pure Yang-Mills theory. Section IV introduces the notion of limited gauge- 
invariance and justifies the gauge-invariance of the results derived. Section V concludes 
with a comparison of the ideas presented here and those found in the literature. The details 
of some of the computations are given in the Appendices 



II. THE DYNAMICS OF THE NEW DEGREES OF FREEDOM 

Consider SU (2) Yang-Mills in 4D Euclidean space-time. Field configurations that satisfy 
equation (1) with d ■ A a = ^f a can be decomposed as: 

Ai{x) = 1 or* + e abc r + rf)(-d,f b + *y , (2) 
(T ' j ~ j ) y 

where f a and are the new scalar (dimensionless) and vector (dimension 1) degrees of 
freedom. These new degrees of freedom satisfy 

W = (3) 

a ^ r^abc , ^abd cd rc ^acd rd rb , ca rd^dbc 

(1+/. f) 2 J J J J 

- r (i + / • f)s bc - r(i + / • fw h ]d,f% = o, (4) 

which make the number of degrees of freedom tally with those of A a . 



Note that the non-linear gauge has a built-in length scale t. In the previous papers, I 
was simply introduced for dimensional reason, i.e., since d ■ A a ~ (length)~ 2 , I is needed 
to make f a dimensionless. Naively, since the entire point of the non-linear regime of (1) is 
confinement, then we should identify £ to be of the order of 1 fermi. But is there a stronger 
basis for t! 

The non-linear gauge condition (see equation (1)) can be written as 

(U 2 5 ab - ge abc A c ■ 8){d ■ A b ) = 0. (5) 



As we will argue in Section V, as we let the coupling constant g runs as given by 



and scans the configuration space, we get to the regime where (d ■ D) ab is singular 
with d ■ A a itself as the zero mode. This happens when g(jj) is strong enough so that 
—g(/j,)e abc A c ■ d is able to lower one of the positive eigenvalues of D 2 to zero eigenvalue. 
Thus, confinement occurs at g(/i ~ Aq C d), and we identify the length scale £ equal to A ^ 

Substituting equation (2) in F" we find that 

F a , v = - g Z%U) + L%,(f; t) + gCfaif; t), (7) 

where Z, L and Q are all functions of f and zeroth order, linear and quadratic in respec- 
tively. Explicitly, the functions are: 

z; u = x abc d,f b d v r, (8) 

L% = R ab (f)(d,t b u - d^) + Y abc (dJ% - d v f%), (9) 

Q% = T abc (fKt c u , (io) 

and the coefficient functions are 

X abc = L_[_( 1 + 2 /. f) e * bc + 25 ab f c -25 ac f b 

(l + f-f) 2 

+ 3e M f d f c - ?>e acd f d f + e bcd f a f% (11) 
R ab (f) = ^^(5 ab + e abc f c + f a f b ), (12) 

yabc = 1 j-=r-[-(f- f)^ + (1 + / • fiff* - (1 - / • f)S°°f b 

(1 + J ■ J) 

+ 3e cad f d f b - 2f a f b f c + e abd f d f c + f a e bcd f d ] ? (13) 

T abc = j T \e abc + (1 + / • f)f b 5 ac - (1 + / • f)f c 5 ab 
(1 + J ■ J) 

_|_ ^abd jd jc _|_ j-a ^bcd jd _|_ ^adc jd y6j (14) 

The derivation of equations (11) to (14) is straightforward but rather tedious. 



Note that in the Coulomb gauge limit (f a = 0), = A ^ and equations (2) to (14) 
consistently reduce to the Coulomb gauge results. 

The pure f dynamics suggests non-perturbative physics because the kinetic term is (<9/) 4 
and its action is ~ p- and infinitely non-linear. This intuition was verified in reference 



(4) where it was shown that a stochastic treatment of a class of classical configurations, 
i.e., those that are spherically symmetric in 4D Euclidean space-time leads to the linear 
potential and the area law behaviour of the Wilson loop. Full quantum treatment of the 
scalar fields, on the other hand, leads to dimensional reduction as the scalar dynamics 
was shown to be equivalent to a 2D 0(1,3) non-linear sigma model via the Parisi-Sourlas 
mechanism. 



Before we proceed in the next section with the discussion on the effective dynamics of 
the "gluon" field (the quotes allude to the fact that we deal with SU(2) and refer to let 
us point out features of the pure f dynamics. 



The class of classical configurations, the spherically symmetric configurations f a (x), not 
only satisfies 

SSf |/ = 0, (15) 



6f'(x) 
but also 

z;m = o, (i6) 

where 

Sf = J d A x\z%Z; v . (17) 

Equations (15,16) follow from the fact that for spherically symmetric fields d^f a = ^^£r 
and X abc is anti-symmetric with respect to b and c. This means that Sf which is positive 
definite, has a broad minimum in the function space of the scalars. It also means that all 
the spherically symmetric scalars are elements of the classical vacuum. It is not apparent 
if the spherically symmetric functions exhaust all classical configurations or the classical 
vacuum. 



What is intriguing about the result of reference (4) is that the linear potential follows 
from a stochastic treatment of vacuum configurations with zero field strength Z° . Is there 
a simple way to understand this result? One argument is that the stochastic vacuum 
configurations reduce the dimension of space-time from 4 to 2 resulting in a linear potential. 
This was confirmed in reference (6) where the full quantum dynamics of f a was shown to 



5 



be equivalent to a 2D 0(1,3) non-linear sigma model. 



III. THE "GLUON" IN THE VACUUM DEFINED BY SCALARS 

In equation (2), we will consider the "gluon" t a in the vacuum defined by the spherically 
symmetric scalars f a (x), where x = (x • x)*. The resulting expansion is not the usual 
background decomposition given by 

A« = A«[f)+t;, (18) 

as can be seen from equation (2) because is coupled in a non-trivial way with f a (x). 
Equation (3) was used in reference (5) to derive the linear potential by assuming a white- 
noise distribution for the spherically symmetric f a (x). Equation (4) on the other hand, will 
lead to the gauge condition 

v; = o, (19) 

which is known in the literature as the radial gauge™. 

Using equations (7) to (14), we find that the effective dynamics of i* in the vacuum 
defined by f a (x) is given by the action 

S e ff{t%f a ) = \J d'x{W\j){d,tl-d^){d,t-d v tl) 

p x x 
+ 2Y ab (f)tpl + 2gV abc (d,tt - dX^K 

+ g 2 T aM (f)t a / v t c /„}, (20) 
where we used equation (19). The coefficients are given by 



K a6 (/~) = R ca (f)R cb (f), (21) 

df d 

dx 

yabffi = Y^(f)Y™\f)^^£, (23) 

U afec (/) = R da (f)T dbc (f), (24) 

rjpafecd rpeab ^j~jrpecd (25) 
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From equation (20), the equation of motion for t a is: 




= w*(f)nX + M ab (fX + Nf u (f, x, d)t b u + P«(/, x) 

+ gU abc (f)d^tl - 9 -W bc {~f){d,tl)tl 
a 2 

2 l v\ l \x l v l v l ix) ~ U ' 



(26) 



where 



M afe (/) 



(27) 



PC(/,x,9) 




6a 



§ a6 )(-)^(x-9) 




(28) 



(29) 



To derive equations (26) to (29), integration by parts and equation (3) were used. The last 
three terms are interaction terms of the vector fields in the presence of the /. The first 
three terms give the propagation of the "gluon" in the vacuum defined by /. The fourth 
is independent and acts like a source term. Appendix A provides the missing steps in 
deriving equations (20) to (29). 

Note the dynamics break translation invariance because of the explicit appearance of x M . 
This is due to the fact that the classical configurations f a (x) implicitly assumes a center, 
the origin. However, this is not a concern because when the averaging over f a (x) is done, 
translation invariance will be recovered. 

At this point, the dynamics seem unwieldy. However, since all spherically symmetric 
f a (x) are vacuum solutions, we should average over all these configurations. As in reference 
(4), we use the white-noise distribution and compute 



where £ was introduced for dimensional reasons. 

The question now is why assume a white-noise distribution for the spherically symmetric 
f a {x)l One reason is that it works, i.e., it yielded the area law behaviour for the Wilson 
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(30) 



loop (see ref. 4). Aside from this, it also means that the "gluon" ft is subjected to a random 
force, albeit, a rather complex combination of white noises. What we have to do is to find 
out if this assumption leads to a testable prediction and the answer is yes and this is the 
subject of the succeeding computations. We will discretize the integral in the exponent by 

/•oo 

/ dsf a (s)r(s) = V f a (s)f a (s)As. (31) 

For each s, the integration in f(s) becomes ordinary integrals involving 

= N- 1 / r 4 e" CTr2 = -(trt)" 1 , (32) 
Jo 2 

poo I 2\m 

j(n,m) = N -l / K ) e -°r* dr (33) 

Jo (l + r 2 ) n 

where m and n are integers, a = ^f- and N = 7i"i<7~^ is the normalization factor, the products 
of which for all s values give the overall normalization M found in equation (30). Equation 
(33) is evaluated using integrals derived from (by differentiating w.r.t. (3 or a) 

roc p—crr 2 _ 

where Reft > and arg o < f (see reference^). The function $ is the error function given 
by 

$(a5/3) = -== / e~ t2 dt, 
V 71 " Jo 

= - \(^f + - ^(^Y + .-.]■ (35) 



The terms to average stochastically are given by equations (24) to (29). There are three 
types depending on the number of derivatives of f a . The terms without derivatives are R, U 
and T given by equations (21), (24), (25), (12) and (14). The results are 

(R ab ) f - = ±5 ab , (36) 
(U afec ) / = (U abc f c ) f = (T abcd ) f = 0. (37) 

The second and third types of terms contain and ^£r^£r respectively. Thus, there is 
a need to define the derivative of a white noise. Since a white noise is continuous but not 
differentiable, i.e., the limit of / ( x +^)-f 0) ag ^ x _^ g ^ oes no ^ ex j s t we need to define 

Ax 
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^£r. This will be done by a suitable smoothening process. Since there is already a natural 
length scale £ in this regime, we will define 

dx ~ £/n ' [6 *> 

where n is a number > 1. Obviously, we cannot take n — > oo, because of the white- noise 
behaviour. 



Using equation (38), all the terms with one derivative, vanish 



<^¥>/=<S->/=0. (39) 



This in turn yields 

(K)f = K)f = °- ( 4 °) 

There are two terms with two derivatives and these are 

<^A/ = (41) 
{Y ah ) } = 5 a \\)^. (42) 
The details of all these computation are shown in Appendix B. 

All these simplify equations (24) and (30) to yield 

q „2 

(° 2 + ^K = 0. (43) 

This means the effective classical dynamics of the gluons show that they are non-interacting 
and have a range of possible masses given by 

3 Th 

At the very least this result verifies the existence of the mass gap for non-Abelian gauge 
theories because the smallest value of the "gluon's" (t^) mass is |(|) = \^qcd- in the 
perturbative regime, the gluon A ^ is massless. 

Is n an integer? Typically, we subdivide a length (in this case the confinement length 
scale t = A q CD into a finite number of equally spaced segments. Experimentally, the way to 



determine n is to compare the predicted mass with the mass spectrum of massive spin one 
"pure glue states". The lowest mass, which gives the mass gap, will determine n m i n while 
the other states will determine the other values of n. Unfortunately, we cannot do this 
comparison because the computation was done for 577(2) and not SU(3) (see Appendix 2). 
But the result will essentially be the same, i.e., m ~ uAq C d, only the proportionality factor 
| will change. 

To conclude this section, we comment on the seemingly paradoxical result of the gluon t a 
acquiring mass and classically losing self-interaction. What was shown is that the effective 
classical dynamics of t a in the stochastic vacuum defined by the spherically symmetric f a (x) 
gives a non-interacting and massive t a . This is not inconsistent with confinement because 
the quantum dynamics of t° in a particular classical background f a (x) necessitates the 
evaluation of the propagator given by 

G bc (x,yJ)=5 ac 5\x-y). (45) 

After computing G(x, y; /), we average over f a (x) to get the stochastic average of the prop- 
agator. This is certainly not equivalent to the propagator of the stochastically averaged 
(M afe D 2 +M afc ) j as can be seen from equation (45). In other words, stochastic averaging does 
not factorize. 

IV. LIMITED GAUGE-INVARIANCE 

At this point, we speculate on the gauge- invariance of the result. The claim is that the 
linear and non-linear regimes of the non-linear gauge satisfy limited gauge invariance. To 
explain what is meant by this, consider first an Abelian theory. Starting from a field A M 
that does not satisfy the Coulomb gauge, one can project to the Coulomb surface to get A M 
or do a gauge transformation to get A' , which satisfies the Coulomb gauge. The results in 
both cases are the same, i.e., 

A; = A, = (V - d^ 2 d v )A v . (46) 
Furthermore, the transverse field configuration A' (and A M ) is also gauge- invariant. 



R ab (f)n 2 + M ab (f) 
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In the non-Abelian case, the situation is more tricky. First, the results of the projection 
(similar to equation (46) only with an additional SU(2) index a) and the gauge transforma- 
tion to the Coulomb surface, (if it exists) are not the same. Second, the projection to the 
Coulomb surface is not gauge- invariant. 

To illustrate limited gauge-invar iance, consider a field configuration such that d-A a = 
e a (x) where e a (x) is infinitesimal and the operator (d • D) positive definite. An infinitesimal 
gauge transformation that will take this field configuration to the Coulomb gauge exists and 
it is 

K = A l~ D f( x ) J d 4 vG bc {x- y; A)(d ■ A%, (47) 
where G ab is the Green function of (d ■ D). Note that if (d ■ D) has L 2 zero mode z a (x; A), 
then the solution given by equation (47) is not valid. Strictly speaking, we need to make 
use of the Green function 

G' ab = G ab - z a (x)z b (y), (48) 

in equation (47) and also impose that J d A xz a {p ■ A a (x)) = 0. This is the reason why we 
require the positive-definiteness of (d ■ D). 

The field configuration given by equation (47) is transverse. It is also gauge-invariant as 
we will show. First 5A'® = D1^(A')a b . The right hand side of equation (47) on the other 
hand can be evaluated by using 5A a ^ = .Djf A b and 

5G ab (x, y\A) = - J d 4 zG ac (x, z; A)8{d ■ D)fG db (z; y; A). (49) 

However, this is too complicated to evaluate. Instead acting by on both sides of the 
variation of equation (47), it is trivial to show that the RHS yield exactly zero^. Thus, we 
have 

(3- D ab (A'))a b = 0. (50) 

Second, since det(d ■ D) is invariant under infinitesimal transformation and since d ■ D(A) 
is positive-definite, then det(d ■ D(A')) is also non-zero and d ■ D(A') is also non-singular. 
Thus, a a = 0, proving that A'^ given by equation (47) is gauge-invariant. Since equation 
(47) is only defined for (d-D) positive-definite, then we were only able to derive a transverse, 
gauge-invariant potential subject to this condition. 
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Now consider field configurations A a (x) such that d-A a = f a (x) and (d-D) ab f b = e a (x) ^ 
0, where e a (x) is an infinitesimal function. Assuming that this field configuration can be 
gauge transformed to the non-linear regime of the non-linear gauge (and as reference (2) 
showed it cannot be gauge transformed to the Coulomb gauge), i.e., A'® = A ^ + D® b (A)p b , 
such that d ■ A' a ^ and (d ■ D ab (A'))(d ■ A' b ) = 0. In this case, 

P a = ~j d A yH ab {x, y] A)[{d-D bc {A)){d-A c )] y , (51) 

where H ab is the Green function of the fourth-order operator 

9 ab = (D ■ d) ac {d ■ D) cb - e acd [d(d ■ A d )} ■ D cb . (52) 

The gauge-potential A'® given by 

A'; = A%- D ab J d*yH bc (x, y;A)[(d- Df d {d ■ A% (53) 

is gauge-invariant. The proof is essentially the same as in the Coulomb gauge only the 
steps are more involved. 

Under gauge transformation, 5A' b = D b ^(A')a c , while 5A b ^ = D b ^(A)A c . Then we act <9 M 
on both sides of the gauge variations of equation (53) and followed by (D ab (A') ■ d). From 
this result, we subtract e abc [d^(d ■ A')]5A'^. The left-hand side is simply 6 ab (A')a b . 

The right-hand side is rather involved because SA^ involves A b while the operator that we 
used involved A'®. Careful evaluation of these terms yield zero^ and a a is then determined 
from 

9 ab (A')a b = 0. (54) 

In reference (12), it was shown in detail that even though the operator d ■ D(A') is singular 
with d ■ A' as its zero mode, the operator 9 (A') is non-singular. Therefore, the only L 2 
solution of equation (54) is the trivial solution and thus proving the gauge-invariance of the 
potential A'* given by equation (53). 

Some details of the above computations are given in Appendix C. 
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We have shown then that the potentials given by equations (47) and (53) are gauge- 
invariant subject to their respective defining conditions. Outside these conditions, these 
potentials are not gauge-invariant in general. Thus, we say that these potentials respect 
limited gauge-invariance. 

The unrestricted gauge-invariance of the Abelian potential defined by equation (46) 
and the restricted gauge-invariance of the non-Abelian potentials defined by equations 
(47) and (53) are reflections of the fact that the gauge conditions used in deriving these 
potentials remove the excess degrees of freedom. These follow from the fact the relevant 
operators - D 2 in the Coulomb gauge of the Abelian theory, (d ■ D) ab in the Coulomb gauge 
of the non-Abelian theory and 9 ab in the non-linear gauge of the non-Abelian theory-are 
non-singular operators. In the case of d 2 it is non-singular without restrictions while 
(d ■ D) ab and 9 ab are only non-singular within limited regimes. In the Abelian theory, if 
is on the Coulomb surface, then + <9 M A will always be off the Coulomb surface because 
□ 2 A 7^ (unless A = 0). In the non-Abelian case, the non- singular ness of (<9 ■ D) ab and 
6 ab in their respective regimes guarantee that the gauge transformed fields in the Coulomb 
gauge and non-linear gauge respectively would no longer satisfy the gauge condition. As a 
consequence, the decomposition given by equation (2), which is only valid if equation (1) is 
satisfied, is not valid for the gauge transformed field. Thus, it does not make sense to ask 
how f a and £° transform under gauge transformation. 

Equation (47) defines transverse vector fields, which can be viewed as the linear regime 
of the non-linear gauge. This is valid when the operator d ■ D ab = n 2 S ab — ge abc d ■ A c is 
positive definite. This happens when the second term is not strong enough to lower the 
lowest eigenstate of n 2 to reach the zero eigenvalue level. And this is true in the per- 
turbative regime where g is rather weak and A ^ is a fluctuation off the trivial vacuum A® = 0. 

As the coupling increases in strength (increasing distance scale), then the lowest mode of 
□ 2 can be lowered down to the zero eigenvalue level. We now hit the Gribov horizon where 
the zero mode of (<9 ■ D) is z a (x). A® still represents fluctuations off the trivial vacuum. 
Here, we cannot define a gauge-invariant potential. There are some claims that non 
perturbative physics is accounted for by this regime. We take a differing opinion because 
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as argued in reference (1), these field configurations do not contribute to the path-integral. 

As we increase further the distance scale, resulting in increasing g, there comes a point 
when (d ■ D) is singular and d ■ A a itself is its zero mode. We now hit the Gribov horizon off 
the Coulomb surface. This is the non-linear regime of the non-linear gauge. In this regime, 
we can define the gauge- invariant potential given by equation (53). Also in this regime, the 
decomposition of the potential in terms of f a and Vt given by equation (2) is valid. The set 
of vacuum configuration is very degenerate, it includes all spherically symmetric f a (x). The 
stochastic dynamics of f a (x) yields the linear potential and generates mass for the gluon t®. 
The full dynamics of f a (x) on the other hand, give the Parisi-Sourlas mechanism. Although 
these are derived from the non-linear regime of the non-linear gauge, we argue that these 
are gauge-invariant results because a gauge- invariant potential (given by equation (53)) can 
be defined in that regime. 

V. CONCLUSION 

In this paper, we argued that the class of spherically symmetric f a (x) (in i? 4 ) form 
vacuum configurations with zero action. Then it was shown that by treating f a (x) as 
stochastic, the effective dynamics of the "gluons" in this background result in the vanishing 
of the "gluon" self-interaction and generation of its mass. 

It is apparent that the class of spherically-symmetric zero field strength (Z® = 0) 
vacuum configurations that produced the linear potential and "gluon" mass is in disagree- 
ment with current ideas on the nature of the QCD vacuum. In the 1970's and early 80's, 
Pagels and Tomboulis^, and SavvidyiS showed that a non-zero field strength configuration 
has lower energy than the zero field strength configuration. Unfortunately, the effective 
Hamiltonian has an imaginary part which signals instability and eventual decay to the 
trivial vacuum. 

About the same time the concept of condensates^, (in particular, the (0 | F 2 | 0) ^ 0), 
supports the idea of a vacuum with non-zero field strength. Obviously, the ideas presented 
here and the condensate result are contradictory. Or are they? 
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Note that the decomposition of the original given in equation (2) shows the scalar 
field f a inextricably linked to the t a and not in the usual background decomposition 
^t, = A® + t a . When we expand the field strength, we get equation (7), which shows 
t a interacting with / in a non-linear manner. The vacuum configurations defined by f a , 
which are the class of spherically symmetric configurations f a (x), have vanishing Z* v . The 
non-zero condensate (0|F 2 |0) involves L a and Q°^ v) i.e., it also reflects the dynamics. 
Thus, the ideas presented here and that of the non- vanishing condensate are not necessarily 
contradictory. 

We also argued the gauge-invariance of the results derived. This is done by showing that 
a gauge-invariant potential, which satisfies the non-linear gauge can be defined within a 
limited regime of potentials. 

Finally, we make the following observation about the non-linear gauge. The non-linear 
gauge condition, as already stated, interpolates between the Coulomb gauge (with weakly 
interacting transverse gluons at the short- distance regime) and the quadratic regime (with 
f a and i° as degrees of freedom, which accounts for confinement and the mass for the 
vector field in the long-distance regime). This assumes that the coupling runs, implying 
that the non-linear gauge is essentially a quantum gauge condition. This is an unusual 
way to fix the gauge. Normally, the gauge condition represents a fixed submanifold in 
configuration space regardless of the distance scale of the physics. The non-linear gauge, 
on the other hand, defines a quadratic submanifold, which, depending on the coupling, will 
highlight either the transverse degrees of freedom of the massless vector field (in the linear 
or Coulomb regime) or the scalar f a and the massive vector field (in the non-linear regime). 
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APPENDIX A: DERIVATION OF THE ACTION AND FIELD EQUATIONS 



From equations (7) to (10), we find 

r ^ rpa rpa 

L-YM - --t^-t^ 

= \ {R ab {d,t - d^)R ac (d,t c u - d v t;) 

+Y abc (dJ% - dJ%)Y ade (dJ% - dj%) 
+2R ab (dX - d v t\)Y acd {dJ c t d v - dj%) 
+2gR ab (dX - d v t\)T ac % 
+2gY abc (dJ% - dJ%)T ad %tl 

+g 2 T abc t b^c T ade t d t e } ^ 

Using dnf b = ^-%r, and equation (19), we find that the Y ■ T term vanishes. Also, the Y ■ Y 
term becomes 

(Y ■ Y) term = 2Y abc Y ade -L-t-t c ^ (A2) 
Substituting these results give the action shown in equations (20) to (25). 

Equation (20) gives the action of the "gluon" t a in a classical background f a (x). In deriv- 
ing the field equation equation given by equations (26) to (29); we made use of integration 
by parts and equations (3) and (19). For example, the variation of the first and second 
terms of equation (20) gives 

first term = R ab D 2 ^ - ±.^^df b 



x dx 
§f c dx x ^ 



5R ab df c A. 

—(- 

5f c dx x' 

1 /In/ h / 1 „f,„ „ n L / 1 



5f c dx y x J M v ' 



second term = ^(§ 6a - § a6 )(-)(x • d)t^ + (-S ba - § a6 )(-)t|l 

Zi X Zi Ou 

+ 2^ { x } g& j 2 5 fd dx*" (A4) 
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APPENDIX B: STOCHASTIC AVERAGES 



We need to evaluate the stochastic averages of the following terms: 

R ab = L^f S °b + fafbs 

(1 + /-/T 



Y ab = 1_. I , 7 f ) 2 < Ml < Ml 6 ab 

(! + /•/) I dx dx 

+ [i + / • / + (/• + (/• /)[i + / • / + (/ • f7\ d £ d £ 

- (i - / • fy cd r^ff b f e T - (/• f)f a f b TT 

ax ax ax ax 

(7 f\ fc^cda fg ^qbe^f df p p. n C cae rb rd df df 1 

-(/•/)/ e f 9 e 9 ^^-(1 -/•/)/ e f f ——\ 

U abc = \ _ (e abc + (2 + f ■ f) e abd f d f c 

(! + /./> 1 

+ (2 + / • f)f a e bcd f d + (2 + / • f)e adc f d f b ) 

T abc d = _ ]_ | [e eab + ^ + ft jag* + fi^ea 

+ e eaf ff f + f e e abf f f + e efb f a f f ][e ecd + (1 + / • f)f c 8 ed 
- (1 + / • f)f d 5 ec + e ec9 f 9 f d + f e e cd9 f 9 + e e9d f 9 f c }} 

S «* = {[(/ • + (l + / • /) 2 / a 5 bc - (i + / • f)5 ab F 

(1 + / ' /) k 

_ (2 + / • /)e aM / d / c - (1 + / • /)/ a / 6 / c + e ac W 

- f a e bcd f d \f- 
dx 

" (! + /•/) 

- 2 (^+r/ 6 )/ c ]^} 
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5S ab df c 

Jflx (1 + f.fy 

+ (1 + / • f~) 3 5 ac 5bd - 2(1 + / • fff a f c 5 bd 

- (1 + / • f) 2 5 ab 5 cd + 4(1 + / • f)5 ab f c f d 

+ (2 + / • f)(l + / • f)e bae f e S cd - (10 + 4/ • f)e bae f c f e f d 

- (1 + / • f)H°*f b f d - (1 + / • fts*/*? 
+ 4(1 + / • f)f a f b f c f d - (1 + / • fiWf'f* 
+ (1 + / • /)e ade / e 5 6c + (1 + / • f)S ac e dbe f e 

- Q e ade f c f e f b -Qe dbe f a f c f e } — — \ (B8) 



-W{[-n/-/)(i+/-/V bd r 



The average of terms without derivatives (B.l), (B.3), (B.4) and (B.6) are of the form 

= AT 1 y (rf/ a (x))F[/]e-Ho°° (B9) 

Using equation (31) and the fact that F[f] is local in x, the averaging will reduce to ordinary 
integration 

(F[f])f = ^J r 2 drdnF(r)e—\ (BIO) 

where r 2 = f(x) ■ f(x), a = and N^ 1 = 7r~i(j + i. At the end of the computation, a — > 
must be taken to get the surviving terms. 

To illustrate the computation, we will show in detail the evaluation of the stochastic average 
of (B.l). 

/ —or 2 
r 2 drcM * + {d 1 * + r a r b ) (Bll) 

Using / dVtr a r b = fr 2 5 ab , 



o ab 



)/ = ~ 2(T+2 ~ i 47r / TT^ dr 







i i 3 j / „ji , / e 



2 



d ^-d e ^ + d (r+^ 
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C O POO 1 /"OO 

" 3 1 (T+T 5 )*} « B12 » 

The first integral is proportional to a~^, the third is a series beginning with cr + 2 and higher 
orders (see equations (34) and (35)) while the second is proportional to <r"t. Thus, in the 
limit a — > 0, we find 

which is the result given by equation (36). The averages of (B.3), (B.4) and (B.6), vanish 
either because (a) they are odd powers in f a or (b) the terms with the highest power of a 
have powers less than — |. For example for (B.4), these are the terms like (1 + / • f) 2 f a f c 5 bd , 
which yield the integral 

°°r 4 (l + r 2 ^ 



(B13) 



o 



Since both the denominator and numerator are of order r 8 , the most divergent part of the 
integral is proportional to 



Jo 



oo 

2 



e~ ar dr ~ a - *. (B14) 



Thus, when we take a — > 0, we get equation (37). 

For terms with one derivative, which are (B.5) and (B.7), we make use of equation (38). 
Now there are two space-time points, at x + ^ and x. Careful evaluation of the terms 
results in the vanishing of the stochastic averages of (B.5) and (B.7). The {f a (x + |))/ 
vanishes while the other term of ^ combines with the other local terms. Fortunately, the 

3 

terms that will contribute a~2 cancel out and the others vanish as a — > 0. These results 
lead to equation (40). 

Finally, terms with two derivatives are (B.2) and (B.8). Amazingly, (B.8) gives zero 
because there is not enough powers of r = (f ■ f) ^ in the numerator to compensate the 
(1 + r 2 ) 4 in the denominator. But the (Y ab )j terms give 
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00 e- ar2 (8 in 16 



+ a+i l (iT^ii rl0 + Y r8 + /j> (B15) 

The dominant terms are the first terms of each integral and the rest all vanish as o — > 0. 
We find: 

3 r? 2 

(Y afc ) = ^(|)(-) (B16) 
APPENDIX C: LIMITED GAUGE-INVARIANCE 

Sincoi^ are unpublished, we provide some of the details in this appendix to prove the 
concept of limited gauge-invariance. 

First, we prove that under condition of the positive-definiteness of (d ■ D), the potential 
given by equation (47) is invariant under infinitesimal gauge transformation. Under gauge 
transformation, the LHS of equation (47) transforms as 

8A 1 ; = Df(A')a h . (CI) 

The transformation of the RHS follows from 

5A; = Df{A)k\ (C2) 

and using equation (49), the RHS gives 

RHS = Df(A)A b + ge M {Df{A)A%J d 4 y G bc (x,y; A)(d ■ A% 

+ Df(A) J d A y J d 4 zG bd (x,z;A)(-g)e def (d-D) f9 A 9 ) z G ec (z,y;A)(d-A c ) y 

- Df(A) J d 4 y G bc (x, y; A) (d ■ D cd A d ) y . (C3) 

Taking the divergence of C.l and C.3, we find that equation C.3 gives zero. Thus, we have 

(<9 • D ab (A'))o b = 0. (C4) 

which is equation (50) of Section 4. Since d ■ D(A) is non-singular, d ■ D(A') is also 
non-singular because det(d ■ D) is invariant under infinitesimal gauge transformation. The 
solution of C.4 then is a a = 0, proving that A'^ given by equation (47) is gauge-invariant. 
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Now consider the potential given by equation (53). Following the same steps as in the 
previous case, we find 

5A'° = Df{A')a b = D ab {A)A b + ge abe {D e /A f ) J d A yH bc (x,y; A)[(d ■ D) cd {d ■ A d ) y ] 

+ D ab J d A yd A zH be (x,z,A)59 ef (z)H fc (z,y;A)[(d- D) cd {d ■ A% 

- Df J d A yH bc (x, y; A)5[(d ■ D) cd (d ■ A%. (C5) 

In the second term in (C.5), we made use of 

6 ab H bc (x, y; A) = 5 ac 5 A {x - y). (C6) 

We will also make use of the following relations in (C.5). 

5[{d ■ D) cd {d ■ A d )} = 6 cd K d , (C7) 

S9 ef = -g e ekl (D l a m A m )d a (d ■ D) kf - g(D ■ d) lk e kfl (d ■ D) lm A m 
- ge ekh [d a (d ■ D hm A m )}D kf 

+ ge ekh [d a (d- A h )]e kfl D l ™A m . (C8) 

Act <9 M on both sides of (C.5) followed by (D a ' a (A') ■ d). When we act this last operation 
on the RHS of (C.5), we take note that A' must be replaced by the expression given by 
equation (53). The LHS becomes (D a ' a (A') ■ d)(d ■ D ab (A'))a b . The RHS is too long and 
rather complicated. However, we can simplify things by neglecting terms proportional to 



[Jd 4 yH(x,y;A)[(d-D)(d-A)] y } 2 A 



because the bracketed term is the infinitesimal gauge transformation p a (see discussions 
before equation (51)) that took A® to the non-linear gauge. The neglected terms are order 
p 2 A and can thus be neglected. With this simplification, we find 

(D a ' a (A') -d)(d- D ab (A'))a b = (D a ' a (A) ■ 8){d ■ D ab (A))A b - g e a ' ac D cd (A) 

x f d 4 yH de (x, y; A) [(d ■ D) ef (d ■ A f )] y d a (d ■ D) ab A b 



+ge abe {D a ' a ■d)(d-D ef A f ) J d 4 yH bc (x, y; A)[(d ■ D) cd {d ■ A% 
+ {D a ' a ■ d)(d ■ D) ab J d A y J d A zH be (x,z;A)59 ef H fc (z,y;A)[(d- D) cd (d ■ A% 

-{D a ' a -d){d-D ab ) J d 4 yH bc (x,y;A)0 cd A d (C9) 
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Now, acting both sides of C.5 by — e a,ac [d^(d ■ A' c ], we find 

- e a ' ac [d^{d ■ A' c )]D ab (A')a b = -e a ' ac [d^{d ■ A c )]D ab (A)A b 
+ e a ' ac [d^d-D) cd J d 4 yH de {x,y; A)[(d ■ D) e f{d ■ A%] 

x D ab A b + g^'^d^d ■ A c )e ahe {D e J V) J d 4 yH bk (x, y; A) [(d ■ D) ki (d ■ A% 

- e a ' ac [d^d-A c )}D ab j d 4 yd 4 zH be (x,z;A) 
x 59 ef (z)H fk (z,y;A)[(d-D) M A e ] 



Adding (C.9) and (CIO), the LHS is 6 a ' b (A')a b . The RHS is zero because: 
(1) The first term of (C.9) and the first term of (CIO) yield 6 a 'A b . 

(2) The fifth term of (C.9) and the fourth term of (CIO) give -9 a ' b A b . 

(3) The fourth term of (C.9) and the third term of (CIO) yield 
56 a ' b (x) J d 4 yH bc (x, y; A) [(d ■ D) cd (d ■ A%. 

Using equation C.8, this term gives 



= {-ge a ' H (D e ™A m )d a (d ■ D) kb - g{D ■ d) a ' k e kbe (d ■ D) hn A m 
- ge a ' ke [d a (d- D) em A m ]D kb 

+ gt a ' kh [d a {d- A h }t kW {D^A m )\ [ d 4 yH bc (x,y; A)[(d ■ D) cd {d ■ A) a ] y - (Cll) 



(4) The second term of (C.9) is cancelled by the third term of (Cll). 

(5) The third term of (C.9) is cancelled by the second term of (Cll). 

(6) The second term of (CIO) is cancelled by the first term of (Cll) 

(7) The third term of (CIO) is cancelled by the fourth term of (Cll). 
Taking all these into account, we find 




(CIO) 




9 a ' b (A')a b = 0. 



(C12) 
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As was shown in^, 9 is a non-singular operator. Thus, the only solution is a a = 0, proving 
that the potential defined by equation (53) is gauge- invariant. 
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